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ABSTRACT 

We show how the existence of a relation between the star formation rate and the 
gas density, i.e. the Kennicutt-Schmidt law, implies a continuous accretion of fresh 
gas from the environment into the discs of spiral galaxies. We present a method to 
derive the gas infall rate in a galaxy disc as a function of time and radius, and we 
apply it to the disc of the Milky Way and 21 galaxies from the THINGS sample. For 
the Milky Way, we found that the ratio between the past and current star formation 
rates is about 2 — 3, averaged over the disc, but it varies substantially with radius. 
In the other disc galaxies there is a clear dependency of this ratio with galaxy stellar 
mass and Hubble type, with more constant star formation histories for small galaxies 
of later type. The gas accretion rate follows very closely the SFR for every galaxy and 
it dominates the evolution of these systems. The Milky Way has formed two thirds 
of its stars after 2=1, whilst the mass of cold gas in the disc has remained fairly 
constant with time. In general, all discs have accreted a significant fraction of their 
gas after 2 = 1. Accretion moves from the inner regions of the disc to the outer parts, 
and as a consequence star formation moves inside-out as well. At z — the peak of 
gas accretion in the Galaxy is at about 6 — 7kpc from the centre. 

Key words: Galaxy: evolution - galaxies: star formation - galaxies: ISM - galaxies: 
evolution 



1 INTRODUCTION 

Star formation is the fundamental process that shapes galax- 
ies into different classes. Although the majority of stars in 
the local Universe are found in spheroidal systems, most of 
the star formation is contributed by disc galaxies of the later 
types (beyond Sb). The key ingredient for star formation is 
cold (star-forming) gas, which is present almost exclusively 
in disc galaxies. However, the amount of cold gas currently 
available in galaxy discs appears rather scant. Kennicutt 
( 1998b I estimated that disc galaxies have current star for- 
mation rates ranging from a few to about ~ 10 M@yr~ l . 
Thus, considering a typical gaseous mass of a few 1O 9 M , 
the gas consumption time scale (i.e. the time needed to ex- 
haust the gas fuel with a constant star formation rate) is 
always of the order of a few Gigayears. This result, known 
as the gas-consumption dilemma ( Kennicutt|1983 1, suggests 



the need for continuous accretion of cold gas onto galaxy 
discs (e.g. San cisi et al.|2008| ). 

Several pieces of evidence show that disc galaxies should 
collect fresh gas from the environment in order to sustain 
their star formation. In the Milky Way, the star formation 
rate (SFR) in the solar neighbourhood appears to have re- 
mained rather constant in the last ~ lOGyr (e.g. Twarog 
1980| |Rocha-Pinto et al.|2000j |Binney et al.|2000[ ) , suggest 
ing a continuous replenishment of the gas supply. Simple 
(closed-box) models of chemical evolution for our Galaxy 
predict too few metal deficient G- dwarf stars than observed 
(|Searle fc S argent ||l972"| |Pagel fc Patchett|[T975} |Haywood| 
2001 1. These observations are easily explained by accounting 
for infall of fresh unpolluted gas (e.g. Chiappini et al.|1997 
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|2001[ > . Observations of Damped Lyman Alpha systems show 
almost no evolution in the neutral gas content of structures 
in the Universe ( Zwaan et al.|2005 Lah et al.|2007 Jorgen- 
|son et al.|2009[ ). Hopkins et al.| ( |2008[ ) pointed out that this 
constancy of gas density can be explained assuming a rate 
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of gas replenishment proportional to the universal SFR den- 



sity. Bauermeister et al. (20101 converged to a similar result 



when comparing the evolution of the molecular gas depletion 
rate with that of the cosmic star formation history (SFH). 
Finally, the derivations of SFHs for galaxies with different 
stellar masses consistently show that late type systems do 
have a rather constant SFR throughout the whole Hubble 
time (e.g. Panter et al. 12007] ). These findings, other than be- 
ing the signature of downsizing in cosmic structures, point 
at a continuous infall of gas onto galaxies of late Hubble 
types. 

The way gas accretion into galaxies takes place is still 
a matter of debate. The classical picture states that ther- 
mal instabilities should cause the cooling of the hot coronae 
that surround disc galaxies and be the source of cold gas 
infall (e.g. |Maller fc Bullock|[2004{ |Kaufmann et al.||2006 |. 
However, recent studies have shown that, due to a combina- 
tion of buoyancy and thermal conduction, hot coronae turn 
out to be remarkably stable and thermal instability does not 



appear to be a viable mechanism for gas accretion (Binney 
et al.|2009 Joung et al.|201l| . On the other hand, accretion 



may take place in the form of cold flows but the importance 
of this process is expected to drop significantly for redshift 
z < 2 ( |Dekel fc BirnboimpjQo] |van de Voort et al.|[20T7 |. 
Observations of local gas accretion at 21-cm emission seem 
to show too little gas around galaxies in the form of HI 
(high- velocity) clouds to justify an efficient feeding of the 



disc star formation (Sancisi et al. 2008 Thom et al. 2008 
Fraternali 2009). However, UV absorption towards quasars 



and halo stars point to the possibility that ionised gas at 
temperatures between a few 10 4 and a few 10 5 K could fill 



the gap between expectations and data (Bland-Hawthorn 



2009} |Collins et al.|2009] |Lehner fc Howk|2011 | 

In this paper, we estimate gas accretion in galactic discs 
indirectly by comparing basic physical properties of galaxies 
today. Our model relies on the existence of a law relating 
the star formation rate density (SFRD) and the gas surface 
density (E gas ) holding at every redshift, i.e. the Kennicutt- 
Schmidt (K-S) law ( |Schmidt||T959"l |Kennicutt||1998b||. This 
approach has similarities with the model of Naab & O striker| 
( 2006[ ) that we describe in detail in Section |4j The paper 
is organized as follows. Section [2] describes our method. In 
Section [3] we apply our model to the Milky Way disc and 
to a sample of external discs. In Section [4] we discuss our 
results. Section [5] sums up. 



2 DESCRIPTION OF THE MODEL 

The evolution of gas and stellar densities (E gas and E») in a 
galactic disc as a function of the lookback time t is described 
by the following equations ( Tinsley|1980| |: 

= +SFRD — Efb (1) 



dt 
dE e 



dt 



= -SFRD + Eft + E 



(2) 



where SFRD is the star formation rate density, Efb is 
the contribution from stellar feedback and E ox t is the in- 
flow/outflow rate from/to the external environment. SFRD 
gives the rate at which stars are formed or, with a sign in- 
version, the rate at which gas is consumed, whilst Efb rep- 
resents the gas density per unit time returned by stellar 



evolution to the gas reservoir in the disc. We treat the stel- 
lar feedback with the instantaneous recycling approximation 
(I.R.A.), meaning that the fractional mass that stars return 
to the ISM at each time is assumed to be constant and equal 
to the return factor 72. Equations |T])-([2| modify as follows: 



dE» 
~dT 
dS g£ 
dt 



= (1 - 7?.) SFRD 

- = -(l-72)SFRD + E c 



(3) 
(4) 



In Section [4. 1| we discuss the effect of a delayed return. 

The parameter 72 in eqs. depends on the IMF 

and the fraction r)(M) of mass that a star of mass M returns 
to the ISM, called the stellar initial versus final mass rela- 
tion. Following the prescription of Kennicutt et al. ( 1994 1 



we assume that high-mass stars (M > 8 Mq) all leave a 1.4 
Mq remnant. Low-mass stars are assumed to leave remnants 
with masses between ~ 0.5 and ~ 1.3 Mg, with a linear de- 



pendency from the initial mass (see Chapter 2 of Matteucci, 
|F.|2001| for details). Assuming a Salpeter IMF, our fiducial 
value for the recycling parameter is 72 ~ 0.30. Assuming 



Kroupa (Kroupa et al. 1993) or Chabrier (Chabrier 20031 



IMFs would make little difference, giving 72 = 0.31 and 
72 = 0.32 (or 72 = 0.46 for a flatter high-mass slope) re- 
spectively With this value the gas consumption timescales 
are extended by a factor (1 — 72.) _1 = 1.4 (see Kennicutt 
et al 1 1994 1. 

The third term on the RHS of eq. (2l is the in- 
flow/outflow term. The combination of eqs. (II I a 
to an expression for this term: 



and leads 



E ex t — 



d£« 
~~dT 



dE e 



dt 



(5) 



In the following we often refer to this term as the accretion 
(or infall) rate. 

Equations |TJ and |2| can be applied to a galactic disc 
assuming that they are valid at each radius, i.e. that the 
disc can be divided into annuli that evolve independently. 
This assumption is equivalent to say that stars that form at 
a certain radius remain at that radius for the lifetime of the 
galactic disc. In Section |4.3| we discuss the effect of stellar 
migration and show that our main results do not change 
significantly. As a consequence of this assumption we can 
integrate eq. |l]) from tf or m (the lookback time at the disc 
formation) to t = to get the current stellar surface density. 

Clearly, in order to solve the above system one needs 
the SFRD as a function of lookback time and radius, i.e. the 
SFH as a function of radius. In the Milky Way, the SFH is 
known with fairly good precision only in the Solar Neigh- 
borhood (e.g. |Rocha-Pinto et al.|2000||Cignoni et al.|2006 |. 
In external galaxies there have been pioneering attempts to 



estimate the SFH as a function of galactic radius (e.g. Gog 



arten et al. 2009 Weisz et al. 20111. In general, however, 



a precise determination of SFH (R) is beyond the capability 
of the current data. Therefore, for the time being, we must 
content ourselves with using simple parametrisations for the 
shape of the SFH. In the following section we describe how 
we estimate the trend of the SFH at each radius in a galactic 
disc. We consider only the disc, leaving aside the bulge/bar, 
which may have formed and evolved differently as it shows 
a markedly different SFH (e.g. Hopkins et al.|2001 l. 
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Figure 1. Evolution of the average SFR, density of the Universe, 
the points are taken from |Hopkins &; Beacom| (I2006J . The thick 
black curve shows a fit with our function f(t) described in the 
text. 



2.1 Reconstructing the SFH of a galaxy disc 

We assume that the SFH in a galaxy disc can be described 
by a dimensionless and positive function / = f(R,t). This 
allows us to write the star formation rate density as: 



SFRD(_R, t) = SFRD(i?,0) x f(R,t) 



(6) 



where SFRD(_R, 0) is the radial distribution of the current 
star formation rate density, which we derive from the data. 
The simplest expression for / one can assume is a first order 
polynomial with time: 

t 



h(R,t) 



h(R) 



(tform >t>0) (7) 



where t is the lookback time and tf OTU1 is the age of the stellar 
disc. With this definition, j(R) is the value of the function 
/ at t = tform and sets the steepness of the SFH. A value 
of j(R) equal to unity represents a constant SFR over the 
disc lifetime, j(R) > (<)1 describes a SFR that increases 
(decreases) with increasing lookback time. The dependence 
of 7 on the radius R allows us to describe different shapes 
of SFH within the disc. In Section [4.21 we discuss the effect 
of taking different forms for /, we anticipate that our main 
results remain unchanged. 

We define the global 7 for a galactic disc as: 



7 : 



2tt J ( f m R 7 (i?) SFRD(_R, 0) dR 
2tt f* m RSFRD(R,0)dR 



(8) 



where the integration is performed out to the maximum ra- 
dius -R m ; the denominator is just the current SFR of the 
disc. This definition allows us to write an expression for the 
global SFR analogous to eqs. (JsJ) and dTb, 



SFR(i) = SFR(0) x 



l + (7-l) 



(9) 



As an illustration of the meaning of our parameter 7 in 
Fig. [T] we fit the average SFR density of the Universe with 
eq. (|9| ). The data points are taken from Hopkins & Beacom 
12006). Here we take tform = 12 Gyr. For the other parame- 



ters we infer SFR(0) = O.O^Mgyr^Mpc" 3 and 7 = 13.4 as 
the best-fit values. This fit shows that the SFR density of the 



Universe as a whole is described by a value of 7 much larger 
than unity, as expected given that the SFR density in the 
past was much higher than now. The shape of the univer- 
sal SFR density is produced by the combined contribution 
of galaxies with very different SFHs. Massive red-sequence 
galaxies dominate the overall SFR at high-z, whilst late- 
type galaxies become relatively more and more important as 
time passes (e.g. Panter et al.|200T Vincoletto et al.|2012 |. 
Although our parametrisation is only suitable for galaxies 
with non-negligible current SFR, taking 7 as the ratio of 
SFR(tf O rm)/SFR(0), clearly red sequence galaxies must have 
7 > 1. 

The above assumptions for the SFRD (eqs. [6] and [7| 
allow us to integrate eq. |3| from tform to to obtain: 

2E„(fl,0) 



y(R) 



(1-11) tformSFRD(7?,0) 



(10) 



where E* (R, 0) is the surface density profile of the stellar disc 
that we derive from the observed surface brightness profile. 
Eq. 10 reveals that our parameter j(R) is ultimately 



a ratio between the stellar density and the current SFRD, 
given a certain formation time for the disc, which is of the 
order of the Hubble time and can be assumed roughly equal 
for all discs. A f(R) close to 1 means that the current SFR 
density extended to an Hubble time produces as many stars 
as observed at that location in the disc. Integrating ■y(R) 
over the disc and making used of eq. Q we find an expres- 
sion for the global 7 as: 



7 ; 



2 M»(0) 



(1 - 11) tf O rmSFR(0) 



(11) 



where M* (0) and SFR(0) are respectively the current stellar 
mass and star formation rate of the entire disc. 

Clearly from eq. one can also describe 7 as pro- 
portional to the ratio between the average past SFR and the 
current SFR. This latter is the inverse of the so-called Scalo 
^-parameter (Scalo 1986) to which 7 is closely related (see 
Section 3.2 1. We also note that any modification of eq. |7| 
with power law terms of ; would produce the same de- 
pendencies for 7 as eq. ( |loj ). We discuss these points further 
in Section l4~2l 



2.2 The star formation law 

The main ingredient of our model is the relation between the 
rate of star formation and the gas density. The empirical star 
formation law is broadly accepted to be a power-law rela- 
tion between the "cold" gas density and the star formation 
rate density ( |Schmidt|1959||Kennicutt|1998b"| ). This relation 
has been tested on global (i.e. averaging gas and star for- 



mation on the whole disc) and local scales ( Kennicutt et al. 
|2007[ ) considering the molecular gas or the sum of atomic 
and molecular gas. Remarkably, it appears to hold across 
various orders of magnitudes, from quiescent discs to star- 
busting galaxies ( |Kennicutt||1998a| |Krumholz et al.]|2011[ ) . 
The dependence on gas metallicity is difficult to investigate 
due to the related dependency of the Xco factor ( jBoissier] 
|et al.|2003| . However, the compilation of azimuthal averages 



of SFRD and E gas presented by Leroy et al. ( 2008 1 for the 



THINGS sample, which includes galaxies with very differ- 
ent masses and presumably different metallicities, seems to 
show that the effect is limited (see Wyder et al.|[2009 and 
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Section 4.4 1. Thus, in the following we assume the following 



star formation law: 



SFRD = AY.I 



(12) 



with N = 1.4 and A = 1.6 x 10 if E Eas is measured in 



M Q pc" 2 and SFRD in A/eyr^kpc" 2 (Kennicutt 1998b I 



Note that the above value of A takes into account the cor- 
rection for the presence of Helium (a factor 1.36 in mass), 
not applied in |Kenni cutt ( 1998b|). 

Deviations from the above relation are observed at col- 
umn densities lower than E gaait h — lOMopc -2 , referred to 
as the density threshold. Below these densities the law steep- 
ens and the scatter increases making it more difficult to de- 



scribe (Schaye 2004 Bigiel et al. 20081. In the context of 



this paper, it may have an effect in the outer parts of the 
discs (see Sections |3.1| and |4.4[) . 



Krumholz et al. (|2012j| have recently pointed out that 



the star formation law expressed in terms of volume density 
is more general than the one expressed in terms of surface 
density. This is very relevant to understand the physical ori- 
gin of the law and to perform numerical simulations. How- 
ever, for our investigation a surface-density law is far more 
practical. Clearly, the choices are equivalent if in the region 
of interest, i.e. within the star forming disc, the scale-height 
of the gaseous disc does not change much as a function of 
radius. The validity of this condition (within a factor 2) is 
well established in our Galaxy, as the disc appears to be flar- 
ing significantly only beyond R ~ 16 kpc ( Kalberla & Dedes 



2008) 



2.3 Inferring the gas accretion rate 

The above assumptions on the SFH and the star formation 
law allow us to infer the local gas accretion rate using eq. 
j5|. If the star formation law is described by eq. (12 I, the 
accretion rate density can be written as: 



E ex t(-R, i) 



(l-1Z)SFRD(R,0)f(R,t) 



SFRD^fl.O) 

NAWttorn 



f(R,ty 



df(R,t) 
dt 



(13) 



where SFRD(_R, 0) is the current total gas density, A and N 
are the parameters of the star formation law and / is the 
functional form of the star formation history. The RHS of 
eq. ( |13[ ) can also be expressed in terms of E gas (7?, 0) (using 
the K-S law) instead of SFRD(i?, 0). This equation is general 
and it could be used if the exact shape of the SFH(i?) were 
known. If / is given by eq. Q, the gas accretion rate density 
becomes: 



t ext (R, t) = (l- 7?.)SFRD(_R, 0) 1 + ( 7 (i?) - 1) 
or, equivalently in terms of the gas density: 



t 

t 



^for 



(14) 



E ext (i?,t) = A(l-^)Ef aa ( J R,0) 



1 + ( 7 (7?) - 1) 



E gas (J?,o) 



MR) 



^forrn 

1-N 



1 + (j(R) 



(15) 



The above equations give the rate at which fresh gas 



must be added to the disc at a certain radius and as a func- 
tion of time to produce the star formation density given by 
SFRD(_R, i). As shown above, the shape of the SFH is regu- 
lated by the value of 'y(R)- For y(R) = 1, SFR and accretion 
are constant in time. More in general, when 7 is close to 1 



the second term on the RHS of eq. (141 is small leading to 
E eKt (R,t) ~ (1 - 7l)SFRD(7?,t). Therefore the gas needed 
is directly proportional to the gas consumed by the star 
formation and this proportionality is simply (1 — TV). For 
7(7?) > 1 the second term in eqs. ( 14 1 and ( 15 1 becomes 



gradually more negative and less accretion is needed. High 
values of "/(R) may result in a negative E cxt (7?, t), hence at 
that time and radius, the amount of gas available in the disc 
exceeds what is needed by star formation. This would be the 
signature of an outflow of gas. 

We define the global accretion rate M ext (t) as the gas 
mass accreted per unit time over the whole galactic disc, ob- 



tained integrating eqs. (141 or (15 I over the disc out to our 
maximum radius 7? m . It is evident that the global accretion 
rate will have dependencies similar to those of the local one 
and the above considerations apply. In particular, for galax- 
ies with 7 of order unity, the global infall rate closely follows 
the SFR at every time M cxt (i) ~ (1 - TV)SFR(t). 

It is worth noting that E cxt (R, t) represents the gas sur- 
face density that has to be added to (or removed from) the 
disc per unit time to produce the reconstructed SFRD(_R, t). 
If we are considering a single annulus, E cxt (7?, i) refers gener- 
ically to the external environment including the adjacent an- 
nuli. Therefore, E ex t(-R, t) coincides with the accretion (in- 
fall) rate at that radius only if there are no radial flows in 
the gaseous disc. This is because, with this method we do 
not trace the location of the gas infall but rather the radius 
at which such gas turns into stars. We discuss this issue 
a bit further later on. On the other hand, M ex t(i) repre- 
sents the total mass needed to be accreted by the galaxy 
disc from outside, at the time t, to sustain the SFR at the 
reconstructed rate. Note however that also in this case the 
definition of outside the disc is merely beyond the maximum 
radius _R m assumed for the disc, which we took to be 5 times 
the scale-length of the stellar disc. 



3 RESULTS 

In this Section we apply the method described above to real 
galaxy discs. We impose the following boundary conditions: 
i) the current stellar density profile E*(_R, 0), ii) the current 
SFR density profile SFRD (R, 0). The latter could be in prin- 
ciple derived from the current gas density profile using the 
star formation law (Section 2.2 1. However here we restrict 



our analysis to cases where the current SFRD is known from 
independent estimates. We first consider the Milky Way and 
then a sample of spiral galaxies with known SFRD profiles. 



3.1 Application to the Milky Way disc 

We assume that the stellar disc of the Milky Way is ex- 



ponential with a scale length hn = 3.2 kpc (Binney & 



Tremainc 



2008) and a maximum extension of R n 



5 h R . 

The normalization of this profile is chosen so that the to- 



tal stellar mass of the disc is 4 x 10 10 M rtDehnen & 
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Figure 2. Current star formation rate density of the Milky Way 
disc estimated using pulsars (open squares), and supernova rem- 
nants (filled squares). These data, originally normalized to the 
Solar neighbourhood value, have been fitted with a function de- 
scribed in the text (solid line) and rescaled so that the integral 
over the disc gives SFR(0)= 3 Afgyr" 1 . The shaded region shows 
the SFRD(i?, 0) one would get from the total (neutral + molecu- 
lar) gas density using the K-S law. The lower and upper bound- 
aries are due to the different determinations of the HI surface 
densities from |Binney fc Mcrrificld (1998) and Kalbc rTa fc Pedes] 
(2008) respectively. 



Binney 19981. The present distribution of the star forma- 



tion rate density with respect to the Solar neighborhood, 
SFRD(_R, O)/SFRD(7? , 0), has been estimated using several 
methods. In Fig. [2] we show the estimates coming from the 
distribution of pulsars (open squares, Lyne et al.|1985 I and 
supernova remnants (filled squares, |Case fc Bhattacharya| 



1998). We fit these points with a function (solid line in Fig. 



2\ of the form: 



SFRD(_R, 0) = SFRD(0, 0) 1 + 



R^ 



(16) 



where a — 3.10 is an exponent defining the tapering towards 
the centre and R* = 1.96 kpc is a scale radius such that the 
peak of the distribution is located at R — R,(a — 1). The 
central SFR density SFRD(0, 0) = 2.3 x lO^Mgyr^kpc" 2 
is normalized in order to have a global current star formar 
tion rate of the Milky Way SFR(0) = 3A'/0yr _1 (see Table 
1 in |Diehl et al.|[2~006[ ). The same normalization is then ap- 
plied to the data points. The shaded region in Fig. [2] shows, 
as a consistency check, the SFRD inferred from the current 
total gas density using the K-S law. The upper and lower 
boundaries refer to two different determinations of the neu- 



tral gas density by Binney & Merrificld ( 1998 1 and Kalberla 
|fc Dedes| ( |2008[ ). The latter is not determined for R < 4 kpc 
The molecular gas density is taken from Nakanishi & So- 



fue ( 2006 ) . On the whole the agreement between the direct 



determinations of SFRD and that inferred by the K-S law 
is very good. In the outer parts the K-S predicts a larger 
SFR than observed, likely due to the low densities and the 
steepening of the law. We quantify this effect in Section 
14.41 Some recent determinations of the current SFR of the 
Milky tend to give values around 1 — 2M0yr _1 (see e.g. 
Murray & Rahman 2010). However, Fig. [5] shows that only 
a SFR(0) ~ 3 M©yr _ is self-consistent with the derivation 
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Radius (kpc) 

Figure 3. Steepness of the SFH as a function of radius for the 
Milky Way disc, parametrised by 'y(R) (see Section |2.1[ l. The 
horizontal line shows the global value obtained for the whole disc. 



of the SFRD(i?, 0) from the gas surface density by inverting 
the K-S law, so we prefer to adopt this value throughout the 
paper. In the following we further fix the age of the stellar 
disc t{ olm = lOGyr or, equivalently, 2f orm = 1.8^] Note that 
taking a larger value for iform, suggested for instance by the 



investigation of Aumer & Binney ( 2009 1 would strengthen 
our results. 

We reconstruct the SFH of the Milky Way disc as de- 
scribed in Section [2.1 1 We find a steepness of 7 = 2.8 for the 
whole disc. The behavior as a function of radius is shown in 
Fig. [3] The parameter ■y(R) appears to be a decreasing func- 
tion of the radius R. Recalling that j(R) is the ratio between 
the initial and the current SFRDs, this trend shows, as ex- 
pected, that star formation in the central regions must have 
been faster than in the outer disc, and it has spread from 
inside out (see also Naab & Ostrikcr 2006). Note that our 
approach is fully axi-symmetric and it does not account for 
the presence of a bar, thus our results for R < 3 kpc should 
be taken with some caution. 

The SFR(t) of the Galaxy reconstructed using the lin- 
ear temporal dependency of eq. |7f is shown in Fig. [4] (top 
panel) (7 = 2.8). We remind that the key parameters to 
obtain this curve are the current SFR and the current 



stellar mass of the disc (eq. 111. Changing these parame- 
ters and the formation time of the disc (iform) would pro- 
duce rather intuitive variations in the value of 7 and the 
steepness of the SFH. Extreme values can be obtained by 
taking SFR(0) = 2M Q yr _1 , M,(0) = 5 x 10 10 M Q and 
iform = 8Gyr on the one hand and SFR(0) = 4M yr _1 , 
M*(0) = 3 x 10 10 Mq and t {olln = 12 Gyr on the other, cor- 
responding to 7 m ax = 7.9 and 7 m i n = 0.8 respectively. 

Having derived •y(R) we can proceed to estimate the 
gas infall rate at each radius and its global value by in- 
tegrating over the whole disc surface using eq. |14| Fig. [4] 
shows the global gas accretion rate necessary to maintain 
the star formation in the Milky Way disc. The accretion 
rate is positive at every redshift starting from 5.4M yr _1 



1 We assume a standard cosmology with f2 D = 0.27 and 
0.73. 
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Figure 4. Top panel: SFR and global gas accretion rate, versus 
time and redshift for the disc of the Milky Way. The dashed line 
shows the evolution of a closed box. Bottom panel: Gas accretion 
efficiency as a function of time, e = 1 corresponds to a complete 
replenishment and thus a constant SFR. 
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Figure 5. Gas accretion into the Milky Way's disc integrated 
over three annuli centered at 3 kpc (short dashed line), 8 kpc 
(solid line) and at 13 kpc (long dashed line) versus time and 
redshift. 



at tfarm and decreasing to 1.3M0yr~ at the current time. 
It follows very closely the shape of the SFR. These results 
are in agreement with earlier determinations based on chem- 
ical evolution models (e.g. Tinsley|1980 Tosi|1988a |. In the 
same figure, we also plot the SFR that our Galaxy would 
have if it evolved like a closed system starting from the same 
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Figure 6. Gas accretion profile at three epochs: lookback time = 
3 Gyr, 1.5 Gyr, and now. Note the very recent development of an 
inner depression and of a peak in the gas accretion distribution, 
now at about i? pea k — 6 — 7 kpc from the Galactic centre. 



initial conditions. A closed system is only able to produce 
~ 8.4 x 10 9 Mq of stars from t = £f orm to t = 0, virtually 
leaving no remnant gas mass at present. On the contrary, 
our model with gas infall predicts a present total gas mass 
of 7.1 x 10 9 Mq, comparable to what observed. The bottom 
panel of Fig. [4] shows the gas accretion efficiency defined as 
the ratio of the accretion rate over the net gas-consumption 
including feedback, e(t) = M e *t(t)/(1 - fc)SFR(t). For the 
whole disc of the Milky Way we obtain £Mw(t) ~ 0.6 at 
t — 0. £mw(£) has a mild dependence on time, reaching a 
value of about 0.9 at z — 1. This shows that the efficiency of 
gas accretion is rather high but not enough to keep the SFR 
constant (e = 1). Thus, gas accretion does not feed the discs 
indefinitely, however it delays the running-out of fuel signif- 
icantly beyond the (closed-box) gas-consumption timescale. 
The value of the current accretion rate and efficiencies are 
rather uncertain as they strongly depend on the assumed 
zero-points. For instance, leaving all the other parameters 
unchanged and taking SFR(0) = 2Mqyt~ 1 would give a 
current accretion of only O.14M0yr _1 . This large difference 
is totally localized to the current time, indeed the efficiency 
of accretion remains always above 0.5 from the disc forma- 
tion to z = 0.1 plunging then to 0.1 at z — 0, see also the 
case of NGC5055 in Section l3~2l 

We now turn to the local infall rate E e xt(.R, t). For vi- 
sualization purposes, we split the galaxy in three represen- 
tative zones: the central part (between — 6 kpc), the solar 
circle (6— 10 kpc) and the outer disc (10— 16 kpc). For each 
region we calculate the local gas accretion rate integrating 
over the annulus surface and we plot the results in Fig. [5] 
We find three different regimes: i) in the central region the 
gas accretion declines very steeply, ii) in the solar circle the 
gas accretion declines by a factor 1.7, iii) in the outer disc 
it increases slightly reaching its maximum at z = 0. The 
behaviour visible in Fig. [5] reflects the simple fact that the 
ratio between the stellar density and gas density decreases 
strongly with radius in the Galaxy. The same kind of trends 
are visible in the SFH reconstructed for the same annuli, 
given that as seen above, E acc (_R, t) ~ (1 — TZ)SFKD (R,t). 

Fig. [6] shows the gas accretion profile as a function of R 
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in the Milky Way's disc now and at two epochs in the recent 
past. The current profile is negative for R < 3kpc, peaks at 
about R = 6kpc and then falls further out. The negative 
value in the central regions shows that there the amount 
of gas is sufficient to sustain the star formation and the 
model would predict outflow (of about 0.2 Afeyr -1 ). This 
inner feature is of recent formation as 3 Gyrs ago the peak 
of accretion was right in the centre. So we conclude that the 
bulk of gas accretion is moving away from the inner disc. 
As a consequence star formation also moves out although 
with a delay with respect to the accretion. Note that the 
shape of the accretion profile at t = f.5Gyr is remarkably 
similar to the current SFRD(i?, 0) (see Fig. GjJ, hinting at 
a delay between accretion and star formation of a Gyr or 
so. Interestingly, this is of the order of the gas depletion 
timescale. 

We stress again that with our approach we are able to 
trace the locations in the disc where the new gas begins to 
form stars. The infall itself could have happened somewhere 
else, realistically more further out and the infalling gas could 
have quiescently flown to inner radii before starting forming 
stars and becoming detectable by our method. 



3.2 External galaxies 

In this Section we extend the application of our method to 
a sample of external galaxies. We consider galaxies from the 
THINGS survey ( |Walter et al.|2008| ) for which radial profile s 
of gas density and SFRD are available ( Leroy et al.||2008 |. 
As for the Milky Way, we assume a common return factor 
1Z — 0.3, an age of the discs of tf OIU1 = 10 Gyr and the linear 
f(t) described in Section 



2.1 We derive the global values 



of 7 both using eq. (Ill and by integrating -y(R) (eq. 10 1 
over the disc. In the former case the values of stellar masses 



and SFRs were taken from Table 4 of Leroy et al. (20081, 



in the latter we used their SFR densities available online. 
The two methods led, in some cases, to significantly differ- 
ent values for 7. This is partially due to the different range 
of integration, Leroy et al. |2008| integrate out to 1.5/225 



while we integrate out to 5 scale-lengths (when the profiles 
extend this far, i.e. for roughly for half of the galaxies) for 
consistency with the Milky Way analysis. 

Fig.[7]shows the derived values for the global 7 as a func- 
tion of the stellar mass (top) and Hubble type (bottom), the 
black star shows the Milky Way. The bars show the ranges 
allowed by the two methods described above. Note that for 
some galaxies also the stellar masses resulting from our inte- 
grations differ from those of Leroy et al.| ( ]2008[ ). The values 
of 7 inferred for these galaxies are all included between 
and 11, with one outlier at 7 = 23, NGC2841, a massive 
Sb galaxy with rather low current SFR (shown only in the 
top plot). There is a clear trend with stellar mass and Hub- 
ble type, showing that big galaxies had much larger SFRs 
in the past. On the other hand, low values of 7, typical for 
smaller galaxies betray a disc continuously forming through 
the acquisition of fresh gas. 

As mentioned, we can relate our 7 to the Scalo b- 
parameter. We use the recent determination of b versus stel- 
lar mass relation from the SDSS ( Brinchmann et al.||2004 l 
and convert it into 7. We plot this relation as a dashed curve 
in the top panel of Fig. [7] The two parametrisations agree 
remarkably well, considering that they are derived in com- 
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Figure 7. The (disc-integrated) parameter 7 (steepness of the 
SFH) for the sample of 21 disc galaxies of |Leroy et al.| l f2008| l as a 
function of the stellar mass (top panel) and galaxy type (bottom 
panel). The star indicates the Milky Way. A common time of 
disc formation of 10 Gyr and the SFH-function f(t) have been 
assumed for all galaxies. The error-bars show the interval allowed 
by two different methods of integration described in the text. The 
dashed curve in the top panel shows the value of 7 derived from 
the Scalo 6-parameter for the SLOAN survey ( Brinc hmann et ah] 
|2004| |. 



pletely different ways. Our values tend to be slightly above 
the curve for masses lower than 10.5 in the log and below 
for larger masses. This is likely due to the shape adopted 
for the SFH, see discussion in Section|4.2| Note that, in this 



respect, Fig. [7] is analogous to Fig. 3 in Kennicutt] ( |1998a[ ). 

For the above disc galaxies we have performed the same 
analysis as for the Milky Way. We derived the radial profiles 
of 'y(R) and £ cxt (-R, i) and integrated the latter to derive 
the global accretion rates. Fig. [8] shows the results for three 
galaxies representative of three classes of stellar masses and 
galaxy types (see labels in Fig.[7|). In general, small late-type 
galaxies tend to have a flatter SFH and an efficiency of gas 
accretion close to 1, see the case of NGC4449 (left panels). 
Galaxies of intermediate types have a shallower SFR and 
efficiency gradually decreasing with time. Finally, massive 
discs with low current star formation rates seem to have ex- 
hausted their capability of acquiring significant amount of 
gas from the environment. They now appear to be consum- 
ing their remnant gas before moving to the red sequence 
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Figure 8. Top panels: SFRs and gas accretion rates as a function of time for three representative galaxies from the sample of |Leroy] 
ct al. (2008). The shaded areas show the uncertainties in the determination of the global 7 described in the text. The open dots are the 
average SFHs estimated by |Panter et a,l.~| ( |2007| > for galaxies of the same stellar masses. Bottom panels: Gas accretion efficiencies as a 
function of time, e = 1 corresponds to a complete replenishment and thus a constant SFR with time. 



(e.g. Cappellari et al.|2011 1. In the same Fig.[8]we also com- 
pare our results with the average SFHs derived using the 
SLOAN survey by |Panter et al.| ( |2007| ) , shown as open cir- 
cles. Panter et al. divide their galaxies in bins of current stel- 
lar masses, those relevant here are 3x 10 10 < M» < lx 10 11 , 
1 x 10 10 < M, < 3 x 10 10 , M» < 1 x 10 10 , respectively ap- 
plicable for NGC 5055, NGC 628, and NGC 4449. The agree- 
ment between our SFRs and their average SFRs is very good, 
except for the galaxies in the intermediate mass class as our 
method tends to find steeper SFHs. Note however that this 
may also depend on the way the stellar masses are esti- 
mated. We used the estimate of Leroy et al. ( 2008 1 who 



simply multiply the 3.6-|J.m surface brightness by a constant 
mass-to-light ratio of 0.5. 

In general, the analysis of these 21 THINGS galaxies 
gives results similar to those found for the Milky Way's 
disc. Gas accretion proceeds almost parallel to the SFR and 
it remains significantly important throughout a large frac- 
tion of the lifetime of a typical disc galaxy. The efficiency 
of accretion seems to decrease fast for large galaxies, while 
small late-type discs remain quite efficient in collecting cold 
star- forming gas from the environment. The current (z = 0) 
values of gas infall that we derive are quite low for most 
galaxies in the sample and a few also have negative val- 
ues showing that the available gas is more than sufficient 
for star formation to proceed. Thus, it would appear that 
these disc galaxies do not need gas accretion to sustain their 
current star formation (see e.g. NGC 5055). This is poten- 
tially a very interesting result but unfortunately, as shown 
for the Milky Way, these values are very uncertain because 
they come from an extrapolation. Moreover, they depend 
on parameters like the stellar masses that have intrinsic un- 
certainties and other parameters that are simply assumed 



(see eq. 14 1. For instance if one increases the formation time 



of the disc or decrease the return factor the value of the 
current accretion rate may increase substantially. What re- 
mains indisputable is the fundamental relation between star 
formation and accretion and the fact that the latter must 
be a major player in the evolution of galaxies from their for- 
mation to the present. There is no way a disc galaxy could 
become as we see it today without substantial gas infall at 
z< 1. 



4 DISCUSSION 

In the previous Sections we have shown how a simple com- 
parison between the stellar density in galactic discs and 
their current SFRD (or gas density) implies the need for 
a large amount of gas infall during their lifetimes. This need 
is rooted in the existence of the Kennicutt-Schmidt law and 
in the assumption that this does not significantly evolve with 
time. Under this assumption we were able to derive the ac- 
cretion rates as a function of galactocentric radius and time 
for a number of disc galaxies including the Milky Way. Here 
we discuss the limitations of our approach and the implica- 
tions of our results. 



4.1 Delayed stellar feedback 

Treating stellar feedback with the I.R.A. saves computing 
time and it allows us to write simple analytic expressions 
for the equations of the model, however the delayed return 
from stars with mass M< SMq may be an important ef- 
fect that needs to be quantified. We added this ingredient 
using the approximation that a star with mass M returns 
a fraction rj(M) of its mass to the interstellar medium at a 
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time tMs(M) (see Kennicutt et al.|1994 1. The stellar lifetime 
(tivrs) are taken from Maeder & Meynet ( 19891: 
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(17) 



given that 4ms (1 Mq) > 10 Gyr, stars below ~ 1 Mq do 
not contribute to the feedback. The difference between the 
I.R.A. and the delayed return in eqs. and @ is in the 
term Eft, that becomes: 

/ioo 
SFRD(t - t MS )»?(M)0(M)dM (18) 

where <f>(M) is the IMF. Given the assumed form of the 
SFRD(J?, t) (eq.[7|, the delayed feedback adds an extra-term 
to eq.[3] Consequently, the expression of y(R) takes the form: 

2 E»(i?,0) 



7df(-R) 



where: 



{1-11+ 2^) tformSFRD(ii, 0) 



1-TZ- 



TDF 



2 T DF 

'form 



t M s(M)r?(M)^(M)dM 



(19) 



(20) 



is a characteristic delayed-feedback time that, with our as- 
sumptions, turns out to be tdf = 3.9 x 10 8 yr (6.0 and 
6.3 x 10 8 for Kroupa and Chabrier IMFs respectively). Eq. 
( 19 1 clearly reverts to eq. ( 10 1 when tdf = 0. Finally, the 



following extra-term: SFRD(i?, 0)(j(R) - l)-^ needs to 
be added to the accretion rate E cxt (R) (eq. \14\ to take into 
account the delay. 

Fig. [9] shows the comparison between the global SFR 
and accretion rate for the Milky Way with I.R.A. and with 
delayed stellar feedback. The values of current SFR and stel- 
lar mass are kept fixed in the two models. As expected the 
global gamma with delayed feedback is lower, 7df = 2.6 (2.5 
for a Chabrier IMF), than that with I.R.A. This because if 
the mass return to the ISM is delayed then less stars need to 
form in order to have the same stellar mass at the end. As a 
consequence also the gas accretion is slightly lower for most 
of the life of the galaxy, reaching however a slightly higher 
value of 1.4Moyr _1 at the current time. Given that the 
extra-term in the accretion formula is positive (for f(R) > 1) 
the efficiency with delayed feedback is always slightly above 
the I.R.A. value. In conclusion, including a delayed feedback 
in our calculations makes little difference and it goes in the 
direction of increasing the efficiency of gas accretion. 

4.2 Parametrisation of the SFH 

The reconstruction of the star formation rate density as a 
function of time requires the parametrisation of the SFH. In 
Section [2. 1| we have chosen a simple first order polynomial. 
We also point out that a more general power-law would give 
similar dependencies for 'y(R). Here we consider three dif- 
ferent parametrisations and show that our results do not 
change significantly. We take the following functional forms: 
i) a piecewise-linear function of time f2(t) that linearly in- 
creases to a maximum SFR at t — t c and remains constant 
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Figure 9. Top panel: SFR and global gas accretion rate, versus 
time and redshift for the disc of the Milky Way with instantaneous 
recycling approximation (thin solid lines) and delayed feedback 
(thick dashed lines). Bottom panel: Gas accretion efficiency as 
a function of time with instantaneous recycling approximation 
(thin solid line) and delayed feedback (thick dashed line). 



at higher lookback times; ii) a triangular function of time 
/s(t) that linearly increases to a maximum SFR at t — t c 
and then decreases to SFR(0) at t — tf orm ; hi) an exponential 
function of this form: /4(f) = e t//r where r is the exponential 
e-folding time of the SFR. With these choices we encompass 
both monotonic and non-monotonic functions. In ii) 7 is 
defined as the ratio between SFR(r c )/SFR(0). 

In Table [l] we show the values of the global 7, and the 
present accretion rate M GX t(0) and accretion efficiency e(0) 
obtained with the four types of f(t) for the Milky Way. We 
took t c = 7.7 Gyr (z c — 1). For the exponential function, 
often used for early type galaxies (e.g. |Bell et aL||2003| ), we 
found a rather long scale-time of r = 8.5 Gyr and a shape 
of the SFH very similar to that obtained with fi(t). Even 
though these three functions depict very different trends for 
the SFH, the disc-integrated 7 varies only by ~ 17% and 
the values of the present accretion rate and e(0) change by 
no more than 13 %. Note that the inverse of the Scalo b- 
parameter would correspond to a piece- wise linear functional 
form /2(t) for t c — ► 0. All the above values are obtained us- 
ing the I.R.A. Using the delayed return with different shapes 
of the SFH can produce different effects. For instance an 
exponential SFH of the type /i(t) maximizes the gas con- 
sumption timescales, which can increase up to a factor 2 — 3 
at the earliest times (see Kennicutt et al.|1994 l. 
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fit) 


7 


M ext (0) 


£ (0) 


linear (1) 


2.8 


1.31 


0.62 


piece-wise linear (2) 


2.5 


1.27 


0.61 


triangular (3) 


2.4 


1.21 


0.58 


exponential (4) 


3.2 


1.49 


0.71 



Table 1. The effect of using different parametrisations for the 
SFH described in the text. 7 is always defined as the ratio between 
the initial over the current SFR except for the triangular function 
where it is the ratio between the SFR at z c = 1 over the current 
value. 
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Table 2. The effect of radial migration of stars on an exponential 
stellar profile and on the disc-integrated 7 for our Galaxy 



4.3 Effect of stellar migration 

In the above we have assumed that radial motions of stars 
in the disc are not significant, here we relax this assumption 



and estimate the effect of a global migration (e.g. Roskar 
et al. 20081 IMinchev et al. 20111. For simplicity we take 



an exponential stellar profile with central surface brightness 
E»(0), a scale length ha, and a maximum extent at radius 
Rm = 5/ir. The stellar mass of the disc is: 



M(£„(0),M = 27rE,(0)ftfl 



1-e" 



Rm 
llR 



(21) 



Consider now a radial motion with a constant velocity v at 
Rm that lasts a time t directed inward or outward. Since 
the initial stellar profile is truncated beyond Rm, all stars 
at this radius will move to the new position R m + vt, which 
represents a new maximum radius and is greater or smaller 
than Rm depending on the sign of v. This condition trans- 
lates to the following relation involving the new (i.e. after 
migration) stellar profile T,' t : 



K(Rm + vt) = T, t (R n 



(22) 



that, together with the assumption that the profile remains 
exponential and that the mass is conserved: 



M'(K(0),ti R )=M(T,*(0),h a ) 



(23) 



allows us to estimate £'*(0) and h' R , the new central surface 
brightness and scale length. 

In Table [2] we report the variation of the exponen- 
tial profile of the Milky Way (S*(0) = 648.2 M pc" 2 and 
flR = 3.2 kpc) in the presence of coherent radial motions. 
When the flow of stars is directed inward, obviously the nor- 
malization of the stellar profile increases and the scale length 
decreases. Eq. |l0| tells us that ~/(R) + 1 oc £*(i?). There- 
fore, the profile of "/(R) will change to 'y'(R) according to 
this relation: 



(24) 



For radii R where £'„ (i?)/£, (R) > 1, eq. (|24| predicts 
an increase in the parameter ^y(R) and vice versa for 
£'»(-R)/£*(7i) < 1. This can be explained considering that a 
radial motion that brings more stars at a certain radius (i.e. 
£„(7?)'/£*(-R) > 1) makes the latter as it had been more 
star- forming. Despite the changes in "/(R) it is remarkable 
that the new global 7 changes very little, see the rightmost 
column of Table [2] showing that our global results are only 
slightly affected by radial migration. 

|Sellwood fc Binney| ( |2002[ ) showed that the dominant 
effect of spiral waves in galaxy discs is to churn the stars 



in a manner that preserves the overall angular momentum 
distribution. Stars move both inwards and outwards without 
any significant radial spreading of the disc or increase in 
non-circular motions. In their Fig. 13 they plotted the final 
distribution of home radii for six different radial bins for our 
Galaxy. The distributions are roughly Gaussian for every 
radial bins, with the largest dispersion a ~ 3 kpc at about 
Rq. So as to give an upper limit on this effect we infer the 
new stellar profile £»(i?)' as the average of that modified 
by inward radial motion with parameter vt — —3 kpc and 
that of an outward motion with parameter vt = 3 kpc. We 
found that the changes in 7(i?) are significant for some radii 
but there is basically no effect on the global value of 7. We 
conclude that both migration and random redistribution of 
stars over the disc do not affect our global results. 

4.4 Comparison to other studies 

It is a matter of debate whether galaxies had their reservoir 
of gas in place from the beginning or they have been har- 
vesting the gas they needed throughout their lives. We find 
that the second must be the case for all late type galaxies. 
A concentration of gas as large as the current disc stellar 
masses would imply an extremely high SFR at earlier times, 
a sudden formation of most of the stars and a steeply declin- 
ing star formation at later times. This is how a red sequence 
galaxy forms and evolves. Galaxies like the Milky Way or 
spirals of later types clearly did not have this kind of star 
formation history. The evidence comes from the reconstruc- 
tion of the star formation histories from the color magnitude 
diagrams (e.g. Harris & Zaritsky 2009[) , from chemical evolu- 
tion models (e.g. Chiappini et al.|2001[ ), and from cosmology 
(e.g. |Hopkins et al.||2008[ ). Our study further supports this 
point making use of the Kennicutt-Schmidt law. 

A possible way out would be to invoke a variation of the 
K-S law with redshift or with metallicity as proposed by e.g. 



Krumholz & Dekel (20111. However, to reconcile the obser- 
vations with negligible accretion at z < 1 the variations in 
the K-S law should be very large. If one required the initial 
gas mass of the Milky Way disc to be as large as the stel- 
lar mass today, the K-S law would imply an initial SFR of 
~ 50Moyr -1 . At that rate, the stellar disc would form in a 
Gyr and the subsequent evolution would be passive (similar 
shape to the closed-box shown in Fig. Therefore to keep 
the current SFR at the observed value the transformation 
of gas into stars should have been an order of magnitude 
less efficient than the K-S law would predict. Moreover, this 
efficiency should have evolved with the galaxy itself in order 
to keep the star formation roughly constant in time. This is 
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Figure 10. Gaseous and stellar mass as a function of the look- 
back time in the disc of the Milky Way. 



in contradiction with the fact that today galaxies with very- 
different masses and thus different disc metallicities roughly 
follow the same K-S law ( |Leroy et al.|2008[ ). Thus, it seems 
unlikely that any realistic evolution of the K-S law can ac- 
count for the evolution of discs without conspicuous gas ac- 
cretion at z < 1. 

Despite the lack of metallicity information in our analy- 
sis, we have tried to quantify the above effect by assuming a 
different star formation law for low density environments. To 
this aim, we broke the star formation law into two parts: we 
assumed that the standard K-S is valid only down to a break 
surface density of 10 Mqpc~ 2 in hydrogen, below which the 
law takes a different form. We took this form from the work 



of Roychowdhury et al. ( 2009 1 who studied a compilation of 



extremely faint (and presumably metal poor) dwarf galaxies 
and found that a normalization of A = 1.0 x 10 -5 and a slope 
of N = 2.47 are suitable for these low density environments. 
This parametrisation makes star formation less efficient at 
low gas densities. We repeated all the above calculations and 
found that the general result is, as expected, that more ac- 
cretion is required, because star formation needs more gas 
to occur at the same rate. However, the differences are not 
big. For instance the amount of gas in the Milky- Way disc 
at all times increases by a factor in the range 9 — 18 % and 
the amount of accretion today only by 5 %. For the other 
galaxies the differences can, in some cases, go up to a factor 
~ 2. The reason for this is that the second term on the r.h.s. 



of eq. ( 14 1 does not change by more than a factor 2.3 below 



the break surface density. 

Another important result of our investigation is that 
gas accretion moves in time from the central to the outer 
parts of the discs. In other words, the inner parts dominate 
the accretion at earlier times, whilst the outer parts domi- 
nate at later times, see Fig. [5] This inside-out evolution is 
also responsible for the production of metallicity gradients 
(see e.g. Matteucci et al. |1999| |Boissier fc Prantzos 1999). 



The progressive shift of accretion may be related to the an- 
gular momentum of the infalling material and it would be 
interesting to check it against the results of ab-initio hydro- 



dynamical simulations (e.g. Piontek fc Steinmetz||20lT and 
references therein). 

Our model also predicts the distributions of stars and 



gas in the disc at all times. We found that the stellar disc of 
the Milky Way remains roughly exponential, at least in the 
inner parts, with a scale-length slightly increasing with time 
from 2.7 to 3.2 kpc. The gaseous disc, initially exponential, 
develops an inner depression with time, as stars are formed 
in the central region and the accretion becomes more and 
more efficient in the outer parts. Fig. |10| shows the evolution 
of the stellar and gaseous masses in the Milky Way disc. The 
cumulative stellar mass increases slowly with time. Roughly 
2/3 of the stars formed at z < 1 as a consequence of gas 
accretion towards the inner disc. On the other hand, the 
gas mass remained almost constant and always rather low. 
Interestingly, the ratio between our initial and current gas 
mass for the Milky Way is the same as that obtained by 
comparing the HI mass in galaxies at z = with that of 
Damped Lya systems (e.g. |Zwaan et al. 2005 Lah et al. 



[2007|. We stress again that the gas total mass in Fig. 10 is 
calculated within i? m = 5/ir, more gas may lay at larger 
radii and flow inwards with time. 



In a recent paper Naab & Ostriker ( 2006 1 built a model 



for the evolution of the Milky Way that has aspects in com- 
mon to our investigation. They consider the growth of a disc 
that goes through an initial merger phase that lasts about 
2 Gyr and then follows a quiescent evolution. This second 
phase is analogous to the phase we are considering in this 
paper. Naab & Ostriker (20061 assume the variation of the 



scalelength of the disc with time as the inverse of the Hub- 
ble parameter, so by a factor about 3 in the last 10 Gyr. 
This assumption together with the K-S law, used to derive 
the SFR, determine the amount of gas accretion onto the 
disc. In contrast to this approach, we derive the SFRD(i?, i) 
from the normalization at redshift zero and use the K-S law 
to derive the amount of gas available at any time and the 
amount of gas needed (gas accretion) to reproduced that 
SFRD(J?, t). Then from this we obtain the evolution of the 
stellar and gaseous discs. 

From the above it is clear that the methods used in the 



present paper and in Naab & Ostriker ( 2006 1 are rather com- 



plementary and it is remarkable that we attain a very similar 
evolutionary picture for the Galaxy disc. The main difference 
is the evolution of the stellar-disc scalelength, which evolves 
qualitatively in the same way but quantitatively much less 
in our treatment, i.e. our disc has initially hn ~ 2.7 kpc as 
opposed to hn ~ lkpc in their model (see their figure 1). 
The SFRs in the Solar neighbourhood are slightly differ- 
ent but both comparable with the determination of |RochaT| 
Pinto et al. (20001. Finally, they also find that the accretion 



rate closely follows the SFR at every time and at present is 
about 2 A/0yr _1 , slightly higher than ours. Our analysis also 
shows similarities to the work of Boissier fc Prantzos| fl999). 
These authors choose a different form of the star formation 
law and derive chemical and spectrophotometric properties 
of the Milky- Way disc. In their case, the accretion-rate is 
fixed and assumed to take an exponential law. Globally, all 
these descriptions give a consistent picture of the evolution 
of a galactic disc, reinforcing the validity of the different 
parametrizations . 



4.5 What type of accretion? 

Our treatment does not include metallicity information but 
it is clear that the infalling gas is bound to be at low metal- 
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licity. The need for metal-deficient accretion has been shown 
and confirmed throughout the last decades starting from pi- 



oneering works such as Tinsley & Larson ( 1978 1 and Tinsley 
( 1980). More recently it has become clear that one needs two 
epochs of gas infall: one at the early epochs and one slower 
and still ongoing, which is the one we are describing here 
( Chiappini et al.|1997 1 . The metallicity of this gas should be 
fairly low, of order 0.1 Zq or less (e.g. Tosi|1 988b). The ori- 
gin of this accretion is still matter of debate and it is beyond 
the scope of this paper. As mentioned, several numerical in- 
vestigations show that cooling of the hot corona may provide 
a reasonable amount of gas accretion onto the Milky- Way 
disc, although at a relative low rate ( Peek et al.|2008 Kauf- 
mann et al.|2009" I . These simulations roughly reproduce also 



the kinematics of this gas ( Kaufmann et aL]|2006 1 and the 
radial distribution of the accretion (Peek 20091. However, 



the whole mechanism of direct cooling via thermal instabil- 
ity has been recently challenged both by analytical analyses 
( Binney et al.||2009| |Nipoti||2010[ ) and by recent grid-based 
simulations (Joung et al. 20111 and the matter is still not 
settled. 

Our accretion is smooth but this is obviously an artifact 
of our assumptions. In principle the main contribution to gas 
infall could be minor mergers with gas rich satellites, which 
would give presumably a more bursty SFH. In this sense, our 
analysis is still applicable if one averages all quantities over 
a time larger than the typical starburst timescale. Current 
HI observations seem to show that the amount of gas that 
can be realistically brought into large galaxies by gaseous 



satellites in the local Universe is rather low (Sancisi et al. 
|2008[ ). However, these estimates are still uncertain and only 
future HI surveys will allow us to determine the accretion 
from mergers with the required accuracy. A recent analysis 
of the contribution of satellites to the gas accretion in the 
Milky Way shows that they would also not reproduce the 
distribution of SFR observed today ( |Peek|2009 1. 

Finally, we found that accretion and star formation pro- 
ceed together. This is clear from eq. |13| and from Figs. [4] and 
[8] This finding may just highlight the obvious point that the 
more a galaxy accretes gas the more stars it forms. However, 
the regularity of this link seems to point at a more funda- 
mental interplay between star formation and accretion (see 
also the discussion in[Hopkins et al.|2008|. In a series of pa- 



pers ( |Fraternali fc Binney 2008; M arinacci et al.|2010|[2011 1 
it has been suggested that gas accretion in discs is funda- 
mentally related to star formation because it is produced 
by supernova (positive) feedback. Supernovae eject (high- 
metallicity) gas clouds out of galactic discs mixing it with 
the surrounding hot corona and causing prompt cooling of 
the latter. This implies an accretion that follows closely the 
star formation in the disc and proceeds from inside out given 
the distribution of cloud orbits ( jFraternali fc Binney] 2006) 



Marasco et al. (20111 applied this model to the Milky Way 
and found a global accretion rate of ~ 2 M0yr _1 . They also 
found that the current accretion profile is an increasing func- 
tion of the radius out to a peak at R ~ 9 kpc followed by a 
sudden drop further out. In this paper, we found a similar 
shape but with a peak located closer in, see Fig. [6] How- 
ever, as mentioned, our approach traces the location where 
the new gas forms stars. So it is conceivable that the infall 
took place at slightly larger radii, and then it flowed quies- 
cently in. It is compelling that two completely independent 



methods as these produce such similar accretion patterns 
and global gas accretion rates. 



5 CONCLUSIONS 

In this paper we have proposed a simple method to derive 
the amount of gas needed for the star formation to pro- 
ceed in a galactic disc, using the Kennicutt-Schmidt law. 
We found that in typical disc galaxies (Sb or later types) 
most star-forming gas is not in place in the disc at the time 
of formation but needs to be slowly acquired from the sur- 
rounding environment. We derived the gas accretion rate 
as a function of time for the Milky Way and other 21 disc 
galaxies from the THINGS sample. We parametrized the 
SFH with a simple polynomial function with a parameter 
(7) that defines the slope of the SFH, 7=1 stands for a flat 
SFH. We summarize our results as follows: 

(i) the disc of the Milky Way as a whole (not only the 
Solar neighborhood) formed stars in the past a rate that 
was ~ 2 — 3 times larger than now; 

(ii) the steepness of the SFH is a function of galaxy mass 
and Hubble type, with late type galaxies having nearly flat 
or inverted SFH; 

(iii) gas accretion is tightly linked to star formation, a 
constant ratio between the two is maintained for a large 
fraction of the life of a galaxy, pointing at a reciprocal in- 
terplay; 

(iv) galaxy discs must have experienced accretion of large 
fraction of their mass (in gaseous form) after z = 1, unless 
the K-S law had a strong evolution with z; 

(v) gas accretion progressively moves from the inner to 
the outer regions of galaxy discs, as a consequence star for- 
mation also proceeds inside-out; 

(vi) the efficiency of gas accretion is less than 1 and 
decreases with time in all galaxies expect some late-type 
dwarfs, as a consequence the gas consumed by star forma- 
tion is not completely replenished and the galaxy eventually 
stops forming stars; 

(vii) the present-time accretion rates we derive are very 
uncertain, but even so there are indications that some large 
spirals may not need ongoing gas accretion. 

An obvious development of this study would be to incor- 
porate metallicity. Our results on the amount of gas infall 
needed for the Galaxy are broadly in agreement with the 
amount of infalling material needed by chemical evolution 
models. Our independent determination of the infall rate as 
a function of time and radius could be taken as an input in 
chemical evolution models to obtain a better understanding 
of the accretion mechanisms that feed the star formation in 
galaxy discs. 
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